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TECHNICAL NOTE 2057 


A METHOD OF COMPUTING SUBSONIC FLOWS 
AROUND GIVEN AIRFOIIS 
By Ate Gelbart and Daniel Reach 


SUMMARY 


An extremely simple computational method la established for 
obtaining circulatory flows around given bodies to a high degree of 
accuracy for flows satisfying the linear pressure and specific-volume 
equation of state. The method depends not on an integral equation but 
on the transformation from the holograph to the physical plane involving 
the determination of an arbitrary analytic function. The dete rminat ion 
of the arbitrary analytic function by elementary means results in a 
close approx im ation of the given body. 


INTRODUCTION 


The real mathematical problem arising from the theory of compressible 
fluids is the examination by the methods of function theory of the 
solutions of the nonlinear first-order simultaneous partial differential 
equations: 


0 = — 
p V 


\ 



• ( 1 ) 


or of the nonlinear second— order partial differential equations 
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where 0 and are the potential and stream functions, respectively, 
and x and y are the independent coordinates in the physical plane 
of the position of a particle of the flov. 

It has "been veil known for half a century that a change of the 
independent variables from the physical coordinates x and y to the 
hodograph coordinates 9 and q transforms the nonlinear system to a 
linear one, that is, equations (l) become 

0 0 = 

0cl “ -T 2 (a)l{r 0 

where 




T 2 (l) = 

p is the density, and M the Mach number. Various methods of function 
theory have been established for the solutions of the linear hodograph 
equations (references 1 to 4) . In particular the theory of E -monogenic 
functions, which is essentially a generalization of the theory of 
analytic functions (references 3 to 5 ), has been used with varied 
success. Analytic functions and E -monogenic functions appear to have 
strong qualitative similarities so that certain solutions of 
equations ( 2 ) "correspond" to well-known analytic functions. 

Because of the difficulties encountered in attempting to solve 
boundary— value problems for E -mono genic functions, "approximate 
q ual itative solutions may be obtained more simply by a method of 
correspondence. The analogous boundary— value problem (the flov past 
a given body) is obtained for incompressible fluids in the hodograph 
plane, and the "corresponding" complex potential in terms of 2 -monogenic 
functions Is then used as a solution to a qualitatively similar problem 
(reference 6 ) . 

This procedure yields a solution In the hodograph plane, which in 
general is not in desirable form. The transformation from the hodograph 
to the physical plane therefore becomes of basic importance. Use is 
made of the transformation of Gelbart (reference 7 ) which was also 
established by Lin and Germain (references 8 and 9) . For the isentroplc 
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• equation of state the transformation has the differential form: 


(3) 

Under the assumption that the pressure and specific-volume relation is 
a linear relation , 1 the holograph equations of the flow, as in the 
incompressible case, are essentially reduced to the Cauchy-Eiemann 
equations, so that 0 and ijr in equation ( 3 ) are essentially harm oni c 
functions. This enables the transitional equation to be -written In a 
relatively simple form: ' ■ 


z = fU> - i f C 7^)- d? (4) 

where £ is the position vector in an auxiliary plane, G( £) is the 
complex potential of the incompressible fluid flow around a circle in 
this plane, with or without circulation, f(£) is an arbitrary analytic 
function of £, and the bar over the integral term indicates the complex 
conjugate. If £ takes on the values at the points of a streamline of 
a flow around the circle in the £-plane, equation (4) is the parametric 
representation of a streamline in the physical plane. The circle as a 
streamline in the £-plane is transformed into the obstacle as a stream- 
line in the z-plane. 

Since for each analytic f ( £) there corresponds a compressible 
flow and vice versa, it is therefore necessary only to determine f ( £) 
such that the circle In the £— plane maps into, a given body in the z-plane 
by means of equation (4) in order to solve the direct problem, that is, 
to obtain the compressible fluid flow around a given body. This paper 
does not attempt to give an exact solution to this problem. The object 
of the paper is to find an analytic, f ( £) by elementary and simple means 
such that equation (4) Is a mapping of the circle onto the given body 
to a high degree of accuracy. 

This is done by considering the following maximum-minimum problem: 
For | £ | = E, to find an analytic f ( £) such that the -maximum value of 


fU) f§~ 4£ -“ (a 



^uch an assumption was first considered by Chaplygin (reference 10) 
and later reconsidered by Karman, Tsien, Bers, Gelbart, and others 
(references 7 and H to 13 ) . 
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shall "be a minimum, where co(£ ) is the analytic fraction that map 3 the 
exterior of the circle onto the exterior of the given tody, the point 
at infinity remaining fixed. The technique of this paper vas less 
specifically set forth in a paper hy Bartnoff and Gelbart (reference 14) 
and consists in expanding each side of the equation 


f(£) /^(o -45 °“ <5) 

in a power series in e i6 and solving for the coefficients of f(£ ) 

hy equat ing the coefficients of the terms with like powers of e i6 . The 
simplicity of this procedure is obvious. 

By sett ing up an integral equation and solving it hy successive 
approximations Bers (reference 15) . is ahle to obtain exact solutions. 
Though the computations are not too long they are many times longer 
than those presented in this paper. 

This investigation was carried out at Syracuse University under 
the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 


SYMBOLS 


0 

If 


x,y ■ 

z 


6 


w 

a 

e 

a 

n» 

1 


potential function 
stream fraction 

Cartesian coordinates in physical plane 
complex variable in physical plane 
complex variable in auxiliary plane 
complex variable in distorted hodograph plane 
magnitude of velocity 

direction of velocity in physical plane , 
velocity vector 

magnitude of distorted velocity 
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<*> 


P 

P 

a 

7 

M 

M 


G(£) 


a 

r 


*($)/»(£) 

«(£> 


E 

6 


A _,3 
n J 1 


0 J V’ 


B n-» D n 


magnitude of velocitj'- at infinity 

velocity of undisturbed incompressible flow 

pressure 

density 

speed of sound 
ratio of specific heats 
Mach number 
stream Mach number 

complex potential of flow around circle in ; plane 

complex potential of compressible fluid flow 
circulation of incompressible flow 
angle of attack of incompressible flow around circle 

angle of attack of compressible flow around airfoil 
analytic functions 

conformal mapping of exterior of circle onto exterior 
of Joukowski profile 

radius of circle in £-plane 

variable between 0 and 2 it 

constant in Joukowski mapping 

auxiliary constants ' , 




T ,T 
1 3 2 


constants 

functions of q. defined by equations (2) 
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as subscripts , particular values of variables 
■with symbol, absolute value 

over term, complex conjugate 
with symbol, total derivative 

THE FUNDAMENTAL RELATIONS 


The four f undame ntal lavs governing tbe steady, two-dimensional, 
irrotational flow of an ideal fluid are: Tbe equation of state, 

p = kp^ v (5) 

where p is .the pressure, P the density, 7 the ratio of the 
specific heats, and k a constant; Bernoulli's equation, 

£_ + P = Constant (6) 

, 2 U P 

the continuity equation, 

nWfJUo (7) 


the circulation equation, 

curl (q) = 0 (®) 

Of these four laws the latter two, the continuity equation and the 
circulation equation, give rise to the potential function 0, given by 

d0 = q cos 9 dx + q sin 9 dy (9) - 

where q is the magnitude of the velocity of a particle at the point x,y 
of the flow and 9 its direction, and to the stream function if, given 

by 


cDjr = -pq sin 9 dx + pq cos 9 dy 


( 10 ) 
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From the fact that d.0 and dijr are exact differentials it follows 

that 


0 = i * 

. P T 


0 = tv 

7 ~ 


1 
P T * 




( 11 ) 


Multiplying each side of 'equation (.10) hy ^ and adding to equation ( 9 ),, 


dz = &- 


19 


(d0 + i di|r) 


( 12 ) 


Equation ( 12 ) is the mapping function from the hodograph plane to the 
physical plane. 

By differentiating equation (12) first with respect to 9 and the 
result with respect to q and equating that with the differentiation 
taken in the reverse order s the hodograph equations of the flow are 
obtained from the real and the imaginary parts : 


<} a = 9. ^ 

y 9 p ¥ q 

0 = — q jLVJLAilr 

1 dq^pty ® 


r 


• > 


( 13 ) 


From equations ( 5 ) and (6) 


the following relations are obtained: 



(14) 


where a is the velocity of sound with 


= *E. 
dp' 



+ 




2 


( 15 ) 
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1 

7-1 


7-1 


P = 


PqI 1 


7-1 qT_ 

2 V 


P 0 U + 


7 - 1 q 


(16) 


M 2_q. 2 __£lp = l 2 

a 2 " P dq - ” a Q 2 - q 2 


( 17 ) 


•where "the subscript; zero indicates the stagnation state. It is 
convenient to no rmali ze the constants so that Pq = &q = 1. This is 

equivalent to introducing the dimensionless variables P jP Q and q_j 
From equation (17), equations (13) become 

*e - p *1 


y'q pq 0 

Upon replacing equation ( 5 ) hy 


J 1 1 > 

S-fx'Hp-pI, 


equations (l4) to ( 17 ) became 


,2 , = _ 


(18) 


( 19 ) 


dp 


( 20 ) 


- p i * 


1 + q 2 - rr - 


( 21 ) 


P 2 = 


1 + q c 


( 22 ) 


M 2 = 


1 + q 


(23) 
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From equations (22) and (23) it follows that 


1 — M c 


= 1 


(24) 


Thus equations (l8) become 


00 = p *q 

d = _ £ •* 

*q q J 


(25) 


There is an elementary transformation that 'will reduce equations (25) 
to the Cauchy-Biemann equations. By a change of the independent 
variables from 0, q to 0,q, where 


1 = 


dq 


(2 6 ) 


^oo q\]l + q c 


and q ro is the magnitude of the velocity at infinity, equations ( 25 ) 
become 






(27) 


After the integration has been performed, equation ( 26 ) becomes 

a -log. (28) 


1 + 


\/i + a 2 


where 


E = 


1 + 


\j! + 1 


2 

00 


q 


00 


(29) 
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From equation (12), 



where Q = 0 + i-^r. From, equations (28) and (22), 

dz = Se lv dO - d 5 ' ' (31) 

2 2JK. 

where w = 9 + ici and 0. is an analytic function of w. Thus, 


-ije- iv dfl (32) 

The function fl is the complex potential of a compressible flow 
in the physical plane or in the holograph plane, depending on in which 
plane it is being considered. However, if ft is considered as a 
function of w, it is analytic and represents the complex potential 
of an incompressible fluid flow in the w-plane. 

Equation (32) will therefore remain invariant under any conformal 
transformation w = w(£). By setting 

- - 1 1o * §§■ (33) 

equation (32) takes on the convenient form 
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’where G(£) = fi[w{£)II and f (£) Is an arbitrary analytic function 
of £ • Since v — Q + iq., the magnitude of the velocity of the com- 
pressible fluid flow can he obtained from, equations- (33) and (28), for 


3 -ie e a = k &L HI 
2 f*(£)' 



( 35 ) 


and 


| G- 1 /f » | 

1 - \& /2f'\ 2 

It is clear from equation (36) that | G* /2f* | < 1. 


( 36 ) 


It appears to he most convenient to fix G(£) to he the complex 
potential of a uniform flow past a circle of radius B ‘ in the £ -plane 
with a constant circulation of magnitude T; that is. 


and 





if n 1 

2« log E 


( 37 ) 


r * ^i^oP Bln a i (38) 

where ctp is the angle of attack and q_ ± is the velocity at infinity 

of the incompressible fluid flow past the circle in the £ -plane. If 
mapping (34) is to keep the point at infinity fixed and is to he such 
that q^ is to he hounded, then f 1 (£) must have the form 

w 0 (39) 

This paper deals only with flows of this type. 

Thus, from equation (35), b Q determines the velocity of the flow 

at infinity. It is easily verified that the proper choice of bp 

insures that the mapping of the circle from the £ -plane to the stream— 
line in the z-plane will also he a closed curve. The other coefficients 
of f 1 (£), and thus f(£), may he so determined that the flow he past 
a preassigned obstacle. 

- It is this procedure that is followed in this paper. 
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TEE DETERMINATION OF THE ARBITRARY ANAI2TIC FUNCTION f(£) 


If the function f(£) has heen fixed in equation ( 3 * 0 ., then each 
point in the plane exterior to the circle E is mapped into a point 
in the z-plane such that- the Telocity of the campressible flow at the 
point z is given by formula ( 36 ) and the boundary of the 
circle ] £ | = E goes into a closed streamline in the z-plane (assuming 
b^ has been suitably chosen) . Therefore* to solve the direct problem* 

that is* to obt ain a uniform flow past a preassigned body* one need only 
determine f(£) such that the transformation in equation ( 3*0 transforms 
the circle ] £ | = E into a preassigned shape, (it is assumed in this 
paper that the confo rmal mapping of .the exterior of a circle onto the 
exterior of any simple closed curve is given. See references 1 6 and 17* ) 

Because the left-hand side of equation (3*1-) is not an analytic 
function* the mapping of the circle onto the preassigned curve is not 
identical* po int by point* with the conformal mapping .of the circle onto 
the preassigned curve. Denote the conformal mapping by co = a)(Q. The 
assumption of this paper* in order to obtain a simple and highly accurate j 
though not exact* means of obtaining f (£)* is that the inverse images 
of the two mappings* the conformal and the one given by equation (3*0 * . 
of the same point on the- given curve subtend a small angle at the center 

of the circle £ = Ee 3 - 5 * 0 = . S < 2 it. This is undoubtedly true for 

airfoil shapes* as has been shown by the computations of this paper. 

That it also is true for curves of circular shapes has been verified in 
the paper by Bartnoff and Gelbart (reference 1*0 . 

On the basis of this assumption the ma x i mu m value of 


- 


1 

IF 



d£ ~ a>({;) 


(**o) 


is pmnn for £ = Be 3 - 5 . By equating the coefficients of like powers 

of e 3 - 5 in z(£) and in cd(£)* and thus determining f(£), it is 
. assumed that the terms neglected in this process will be small. This 
has been verified for the cases chosen in this paper and in reference 1 **. 

The process here mentioned will now be applied to the case of a 
symmetric Joukowski airfoil for a given angle of attack ct and a flow 
whose speed at infinity .is q^. The fact that a symmetric Joukowski 

airfoil is chosen instead of a general airfoil is by no means a necessary 
restriction. The same method with the same simplicity can be applied 
to all airfoil shapes. 
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To' satisfy the Kutta-J oukowski condition at the trailing edge, the 
flow past the circle £ = Be i5 with an angle of attack 0 ^ or a 
circulation 


• r 4«q. B sin a. 

-L,oa 1 


(41) 


is given by 




■ia^ J 2ir 


IT _ te 

_ _ iog 


-icu 


B 


(te) 


The stagnation points on the circle occur at 6 = 0° and 6 = l8o° + 2a., 
Taking the derivative of G(£) and substituting equation (4l) for f. 


G' 


<£> = a i,» e 


2iB sin a* 


^ e “ ict l ^2 e “ 21a i, 


(43) 


Squaring each side of equation (43), 


t>' (£ 0 2 - *1, > 21a f 0 +A lf + A 2^ + A 3 ^ + A h^j 


where 


A q 1 


A x = 4iBe iai sin 
0 2ia. 

A 2 = H3ETe 1 (l + 2 sin 2 a^) 
A 3 = — 4iB^e^^ a ^sin 


, 1 , 4ia., 

\ = E^e 1 


(te) 
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Eecall that f'(0 = Let 

n=0 = 


i -V hi 

f'(£) »n 


11=0 £ 


Hence 


so that 


\n=0 i ) \n=0 £ 


1 = HqBq 

0 “ t^o +. • * * + t 0 B nJ n > 0 


Equations (47) can he solved for B^: 

B 0 =i- 
■ t 0 




B n = “htt^O + * * * + b l B - 


'n-l) 


s' 


Por later use B-^ and B 2 in terms of the nhkno'wn constants h^ 


given: 


B X = 


b 2 = 


D 0 

h 2 _ *>1 


v 


+ 

-t 2 -u 3 

B 0 B 0 


J 


Now it is possible to write 


fcii SL = a -1 

■f’(5) 1^0 n t n 


(45) 


(46) 


(47) 


(48) 

are 


(49) 


( 50 ) 
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where 


D n = AcPn + • ' * + A 4 B n-4 


(51) 


Thus, 


P l&( 

U f 


m 


( 5 ) 


d£ = 


2^ 21a i / ' 


D. 


' 11=1 n£ 


and 


(52) 


f(0 = + b Q £ + h log g -X — 

n=1 nf 

On the circle £ = Re 18 ,. equation (34) hecames 


(53) 


z = h_-L + hoRe 18 + h-L log R + h 1 (i5) e -in8 - 


n=l 


hR n 


°° D. 


2 e 210 ^ / ’ 

W 0-15 + 5 X lo S H " »i(») -H ^ 

4 ' • n=l nR n 


n+1 Q ln6| 


\ 


(54) 


• The Joukowski transformation that maps the exterior of a circle of 
radius R onto the exterior of an airfoil is given by 


0) 


(5) + 

£ - 6 


= -€ + £ + Z1 

n=l 


(R ~ 6) 6 

£ n 


2 n— 1 


(55) 


where e is the center of the circle that is mapped onto a strip. On 
the circle £ = Re 18 , equation (55) "becomes 


>2 n— 1 


CD 


( E e«) = -€ + Eei s + f: fe ~ *> £ e — InS 

n=l B n ' 


(56) 
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Equating coefficients of like terms in equations (5^) and ( 56 ), the 
following simultaneous equations for b^ are obtained: 


p 

q.± M e _ 

-e = b_-^ + b^ log B Z—r I°S ® 


(57) 


E = Bb r 


0 = bi + 




2 210^ 


5-i,« 


■j Dp 

lffi 2 

t 

i 


2 2ia 




(58) 


(59) 


ft =jflf - -I 2 - 

E 


E 


2 2icL« 

a 1 • e 

-±i2 ED, 


( 60 ) 


and 


(B - e) 2 e n 1 _ b n+l n > 2 


E n nR : 


n 


( 61 ) 


Since D n is a function of b Q , b x , 
(59), and ( 60 ) can be solved for 1 )q, 


b n only, equations ( 58 ), 


b„. and bp. Once these are 


obtained, equation (57) can be trivially solved for b_ r Equation (61) 
can also be trivially solved for b n , n = >3. It remains then only to 
solve equations ( 58 ), ( 59 ), and ( 60 ) for b 0 , b x , and b 2 . 

Equation (59), it should be pointed out, is the condition on’ b x 
that the rasping of the circle £ = Ee i5 hy emotion OW *= a closed 

Eguations ( 58 ), (59), ^ (60) ore definitely not linear emotions 
“ . But if, from, equation ( 35 ), the value of q i ^ 00 ln 
°is substituted into the three equations they become linear 

_ — - * 1.1 L 


curve. 


in b Q or b 
terms of 


in P 


0> 


and bp, -where P n is the absolute value of b Q , 


that 


is * P 0 = Po 
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For the point at infinity equation (35) gives rise to the relation 


e* a = 


. — ia, 
q. e ^ 


2 Dr 


(62) 


or 


g i ( CCr-CXj) i ( Ctr—Qj, \ 

= - -h « v 1; = 2Fb.e V V 


A,« = f V 


(63) 


where a is the angle of attack of the compressible fluid flow past the 
airfoil, and 


p = I = 
K 


1+ { 1 + k 


( 6 k) 


0 < P < 1 


(65) 


Since ^ and P are real, the argument of h 0 must he a - clj, 
i^ar-a^ 


or h Q = P Q e 


. Equation (63) then becomes 
li,. * 2^0 


( 66 ) 


After substituting equation (66) and Dq = ~ into equation ( 60 ), 

D« 


b 2 = _(E - e) 2 - P 2 p 0 R 2 e i ^ CW1,1 ^ 


( 67 ) 


Now 


D 1 “ -^l + A 1 B 0 

_ ^ _ ^-iBe ^ sin ccj 

So 2 


( 68 ) 
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The substitution of this and equation (66) into equation (59) yields 

bj — P 2 e 2ia 'b’i — 4iP^e ia Po sin = 0 (69) 

The complex conjugate of equation (69) results in 


^ = P 2 e _2iCL b 1 - IdP^e'^Q sin a ± (7°) 

and by substitut ing equation (70) into equation ( 69) equation (69) 
becomes 

b - p\ 1 + 4iP i, Ee ia P 0 sin - l4-iP 2 Re 1 ^ 0 sin a ± = 0 (71) 

Hence 


4iP^e la sin a. 


*1" 


(72) 


1 + P 


Proceeding as before, 

D 2 = AqB 2 + AjB-l + 


_ - t- 2 

+ ^ + 


4iRe 


— iou, 


sin 


\ 2 ”o 3 


V 2 


2R 2 e _2iCti (l + 2 sin 2 ^) 


(73) 
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The substitution of this as well as equations (72) and (67) into 
equation (58) yields a linear equation in 0^: 

- p/O^E 2 + P 2 e 2ta (E - 0 2 + - 

l6P^2e 1 ( a ^i)p 0 

sin 2 ^ 

(1 ♦ p 2 ) 2 . 

l6P 1< R 2 e i ( CCr<Ci )pQ sin 2 ^ 


2P 2 R 2 P o e 1 ^ ar ^ 1 ^l + 2 sin 2 ^) 
Thus the solution for 0 Q can be written 


( 7*0 



K 2 _ p 2 e 21a 


ft ~ e) ! 


h P 2 ^ 2 

/ 2 P sin 

2 

i 1 - p ) 

1 \ 1 + P 2 y 

. 


(75) 


In determining the coefficients of f(£) it is desirable to do 
so in terms of q w and a, which are to be considered as preassigned. 

The quantity must then be considered as an unknown. The four 

equations for 0 Q , b.^ b^ and a are then equations (58), (59), 

(60)_, and (.63). Since 0 Q must be real, Oj can be obt ain ed by 

setting the imaginary part of equation (75) equal to zero: 
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or 


-E 2 sin — a,^ — P^CE — e) 2 sin + a^= 0 


(77) 


Hence 


sin (a - c^) = p 2 ( R _ e)g = 
sin 


— Q 


(78) 


It is seen from equation (78) that 0 is positive*. Prom, equation (78) 


sin cos a, — cos sin a, 

sin a„ cos a, + cos a.^ sin a 


- 0 


(79) 


and 


tan 


tan 


-SL- = q/i + 

a i \ ^ aa a i ) 


■ (80) 


Then for a,. 


tan a. = + ^ tan a 

1 1-Q, 


(81) 


where 


Q . , 


(82) 


it 


Since 0 = a < always has a solution between 0 and g 
and cc^ ^ a. When a = 0, = 0; 

To obt ain b^ n ^ 3, equation (6l) is solved: 

T. n = -<n - X)(E - 


a = 3 


(83) 
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Final l y, 


*■(£) 5a 

3 Z* 




\ 00 

l n + — + — - 

0 i *2 


( E - e )2( n - l)€ n-2 


r n= 3 


5 


n 


b, 1, 


*0 + f + "I + (E " e)i 
5 6 


£ 2 ( C -*) 2 


(84) 


Also, 


f(5) = ^ ^ log ^ - — (B — lL , (g. - g ) 2 


c 


i -c 


(85) 


where b_^, b^, b^, and bg are given "by the equations (57), (75), 

(72), and (67). 

The fact that a J oukowski airfoil was chosen is in no way a limita- 
tion of the method. The same procedure could he followed with the same 

OO . 

simplicity for the conformal mapping function «(£) =a_^ +21 A n /r 

n=0 

of the exterior of a circle onto the exterior of an arbitrary airfoil 
shape, the point at infinity remaining fixed and the value of the 
derivative at infinity being bounded. 

If nonsymmetric Joukowski airfoils are considered, there is no, 
essential change in the representation' of b . Here, e would be complex, 

but at no point in the work has it been essential that g be real. 

It should be remarked that the method of this paper does not yield 
an exact solution of the direct problem, but for practical purposes' 
gives -a very good approximation to the exact solution. However, the 
results obtained by this method are exact for the flow aro und a .body 
that can be computed readily and for which there is a. precise mathe- 
matical expression, such as the Joukowski airfoil used in the example. 

Bers (reference 15) approaches this problem by means of an integral 
equation, which yields exact results whenever the iteration process is 
valid, but the computations are much more involved. 
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ESSUITS OF COMPUTATIONS 


Computations hare "been carried out for a given Joukowski profile 
having the following geometric characteristics: e = 0.15 and R = 1.15* 

The constants of the transformation for free-stream Mach numbers of 0.5 
and O.685 are given in table 1 . The velocity distributions calculated 
for the actual profiles for various angles of attack and the same two 
values of the Mach number are given in tables 2 and 3 and plotted in 
figures 1 and 2 . Tables 4 and 5 present the pressure distributions 
corresponding to the velocity distributions of tables 2 and 3 • 

The velocity computed at the trailing edge (£ = B) is zero unavoid- 
ably because G' (£) is zero at this point in the formula, 


1 = 


G '(0 

f* ( 5 ) 

s T (S) 

2 f'U) 


2 


( 86 ) 


and f (£) is not zero. The exact f(£) which maps a circle onto a 
Joukowski profile would have a zero derivative, but the one obtained by 
the present method is only an approximation. Accordingly, the data 
for' 5 = 0 ° and 360° have been omitted from tables 2 to 5 and 
figures 1 and 2. 


ACCURACY" 


The methods of this paper are primarily for computational purposes 
and are to be Judged by the simplicity and accuracy of the computations . 

The body for which ~the exact flow is. given by the work of this 
paper is represented by the equation 


z = co(£) 



(87) 


where cd(£) is the mapping .function of the given body. The accuracy 
then depends on how smal 1 the second term of the right-hand side of 
equation (87) is on the circle £ = Ee i& . The ordinates for the given 
body and for the one obtained from equation (87) are given in table 6 
for M ra = O.685 and a. = 0 °; the corresponding profiles are given 
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in figure 3* For this case, the series term in equation (87) is an 
alternating series of decreasing terms which begins with the term. 

• involving . The error in neglecting the remainder is then always less 

than its first term. The maximum value of the term involving D did not 


exceed 0.0002. 


The sum 


11 D , 
v — n+1 

ii=2. nl n 


on 


t = Ee i& 


was computed and nowhere 


exceeded 0.022. 


The mapping function f(£) to "be determined for a given flow around 
a given body not only depends on the given body but also on the free— 
stream velocity and the angle of attack. This implies that the body 
obtained from mapping (87) varies slightly for different free-etream 
flow, although <u(£) remains fixed. 

In table 3 and figure 2 of the present paper and also in the paper 
by Bers (reference 15) are data for the velocity distribution of the 
flow around a Joukowski airfoil,, for which e = 0.15, E = 1 . 15 , 

Mo, = O.685, and = 2° 27* • Although a precise comparison cannot be 

made since there is no precise way of fixing the points of comparison, 
the two results appear to compare quite favorably . 


Syracuse University 

Syracuse, E. Y., October 27, 19^7 
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M„ = O.53 q. = 0.57743 P = 0.2679; 6 = 0.15; E = 1.15 



0° . 

5° 

H 

O 

O 

15° 

20° 

a 

p o 

0° 

1.0976 

40. 29* 
1.0934 

8° 59' 
1.1063 

13° 31* 
1.1243 

18° 5* 
1.1431 

,kl/ E 

0.0000 

0.0255 

0.0515 

0.0780 

0.1048 

|l + i 2 | /e 2 

0.0788 

O.O785 

0.0794 

0.0807 

0.0821 

^•1,00 

0.5881 

O.5858 

0.5928 

0.6024 

0.6125 

M„ = 0.6853 < 

3 = 0.94023 P = 

co 

0.3963; « 

’ = 0.15; 

E = 1.15 


0° 

ro 

0 

ro 
— ] 

5° 

10° 

15° 

a 

p o 

0° 

1.2403 

2° 43* 
1.2390 

' 3 0 56’ 

1.2351 

7 0 54* 
1.2619 

11° 55* 
1.2966 

N/ E 

0.0000 

0.6405 

0.0585 

O.1190 

0.1823 

I 1 + \\I S2 

0.1949 

0.1946 

0.1940 

0.1982 

0.2037 


HjOS 


0.9831 


0.9820 


0.9789 


1.0002 


1.0277 
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TABLE 2 .- VELOCITY DISTRIBUTION 


DISTRIBUTION qyq.^ FOB VARIOUS ANGLES OF ATTACK 
[m^O.5; 6 = 0 . 15 ; E = 1.15] 



0° 

5 ° 

10° 

15 ° 

0.875 

0.871 

0.867 

0.855 

.889 

.895 

.901 

.902 

.912 

.926 

.941 

• 950 

.943 

.967 

• 991 

1.009 

.982. . 

1.016 

1.050 

1.079 

1.027 

1.073 

1.119 

I.161 

1.077 

1.137 

1.199 

1.255 

1.129 

. 1.207 

1.287 

1.362 

1.182 

1.283 

1.385 

1.484 

1.234 

1.362 

1.492 

1.621 

1.281 

1.443 

1.611 

1.777 

1.319 

1.526 

1.741 

1.960 

1.342 

1.608 

1.887 

2.177 

1.339 

1.684 

2.054 

2.449 

1.285 

1.744 

2.249 

2.809 

1.126 

1.749 

2.458 

3.297 

.734 

1.550 

2.542 

3.839 

.000 

. 875 ' 

1.955- 

3.526 

.734 

.000 

.792 

1.813 

1.126 

.560 

.000 

.609 

1.285 

.852 

.431 

.000 

1.339 

1.002 

.674 

.343 

1:342 

1.077 

.816 

• 553 

1.319 

1.107 

.98 9 

' .685 

1.281 

1.110 

.942 

.768 

1.234 

1.096 

.960 

.817 

1.182 

1.071 

.961 

.844 

1.129 

1.040 

.950 

.855 

1.077 

1.005 

.934 

.856 

1.027 

.970 

.914 

.851 • 

.982 

.938 

.894 . 

.844 ' 

.943 

.910 

.876 

• 837 

.912 

.888 

.864 

‘.834 

.889 

-.873 

• .858 

.837 

.875 

.869 

.864 

.850 



1.195 

1.328 

1.^30 

1.575 

1.743 

1.941 

2.179 

2.478 

2.874 

3.444 

4.332 

.5.685 

6.218 

3.423 

1.362 
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TABLE 3.- -VELOCITY D3STRIBUTI01I <3./^ EQR VARIOUS AEGEES OF ATTACK 
[M w = O.685; € = 0.15; B = 1.15] 


'\ a i 

6 X. 

0° 

t— 

CM 

0 

CM 

5° 

10° 

15° 

10 

O.876 

O.87O 

0.866 

0.858 

0.843 

20 

.894 

.897 

.896 

.904 . 

.907 

30 

.919 

.929 

• 931 

•9^9 

.963 

4o 

• 953 

.970 

.976 

1.004 

1.028 

50 

.996 

1*022 

1.029 

1.069 

1.106 

60 

1.045 

1.079 

1.092 

1.146 

1.197 

70 

1.100 

1.145 

1.162 

1.234 

1.303 

80 

1.159 

1.218 

1.241 

1.334 

1.425 

90 

1.219 

1.295 

1.326 

1.445 

1.566 

100 

1.279 

1.375 

l.4l6 

1.569 

1.727 

no 

1.333 

1.456 

1.510 

1.707 

1.913 

120 

1.377 

1.535 

1.605 

1.860 

2.133 

130 

1.403 

1.608 

1.700 

2.034 

2.401 

140 

1.400 

1.666 

1.789 

2.235 

2.744 

150 

1.338 

1.691 

1.860 

2.475 

3.212 

160 

1.157 

1.627 

1.864 

2.746 

3.902 

170 

.731 

1.320 

1.637 

2.908 

4.875 

180 

.000 

.583 

.893 

2.280 

5.092 

190 . 

.731 

.224 

.000 

.856 

2.453 

200 

1.157' 

.741 

.565 • 

.000 

.684 

210 

1.338 

1.013 

.874, 

.442 

.000 

220 

1.400 

1.146 

1.036 

.696 

.358 

230 

1.403 

1.205 

i.n6 

.846 

.576 

240 

1.377 

. 1.219 

i.i48 

.931 

.714 

250 

1.333 

1.207 

1.149 

• 977 

.800 

260 

1.279 

1.179 

1.132 

.99^ 

.851 

270 

1.219 

1.140 

1.102 

.992 

.877 

280 

1.159 

1.096 

1.066 

.979 

.887 

290 

1.100 

1.051 

1.026 

•959 

.885 

300 

1.045 

1.007 

.987 

.936 

00 

CO 

• 

310 

* .996 

.967 

.952 

.913 

.868 

320 

.953 

.932 

.920 

.893 

.860 

330 

.919 

.905 

.896 

.879 

.856 

340 

.894 

.885 

.881 

.873 

.859 

350 

.876 

.876 

.877 

.878 

.873 

■ 
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TABLE 4— PRESSURE DISTRIBUTION 


*/*0 


FOR YARIOUS ANGLOS OF ATTACK 


.5; e = 0.15; R = 


\ a i 

8 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 . 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 


0 ° 


1.120 

1.124 

1.130 

1.139 

1.150 

1.163 

1.178 

1.194 

1.211 

1.228 

1.244 
1.257 
1.265 
1.264 

1.245 

1.193 

1.086 

1.000 

1.086 

1.193 

1.245 

1.264 

1.265 
1.257 

1.244 

1.228 

1.211 

1.194 
1.178 
1.163 

.1.150 

1.139 

1.130 

1.124 

1.120 


5° 


1.119 
1.126 
1.134 

1.145 
1.159 

1.176 

1.196 

1.219 

1.244 

1.272 

1.302 

1.333 

1.365 

1.395 

1.419 

1.421 

1.342 

1.120 
1.000 
1.051 

1.115 

1.155 
1.178' 

1.187 

1.188 

1.184 

1.176 

1.166 

1.156 

1.146 

1.137 

1.130 

1.124 

1.120 

1.119 


10 ° 


1.118 

1.127. 

1.138 

1.152 

1.169 

1.191 

1.216 

1.246 

1.281 

1.320 

1.366 

l.4i8 

1.479 

1.551 

1.639 

1.736 

1.776 

1.508 

1.100 

1.000 

1.030 

1.073 

1.106 

1.126 

1.138 

1.143 

1.144 
l.l4l 
1.136 
1.131 

1.125 

1.121 

1.118 

1.116 

1.118 


15° 


1.115 

1.127 

i.l4l 

1.157 

1.178 

1.204 

1.235 

1.272 

1.317 

1.370 

1.433 

1.510 

1.606 

1.732 

1.905 

2.151 

2.432 

2.268 

1.448 

1.060 

1.000 

1.019 

1.050 

1.075 

1.094 

1.106 

1.112 

1.115 

1.116 
i.n4 

1.112 

1.111 

1.110 

1.111 
i.n4 


20 ° 


1.111 

1.126 

1.142 

1.161 

1.185 

1.215 
1.260 
1.297 
1.352 
1.419 

1.502 

1.607 

1.746 

1.938 

2.226 

2.69V 

3.432 

3.727 

2.215 
1.272 

1.036 

1.000 

1.013 

1.036 

1.056 

1.072 

1.083 

1.090 

1.095 

1.097 

1.098 

1.099 
1.101 
1.104 
1.109 












w 


1.402 

1.439 

1.479 

1.521 

1.564 

1.603 

1.636 

1.656 

1.653 

1.607 

1.477 

1.213 

1.000 

1.213 

1.477 

1.607 

1.653 

1.656 

1.636 

1.603 

1.564 

1.521 

1.479 

1.439 

1.402 

1.370 

1.343 

1.322 

1.306 

1.295 


2° 27' 

5° 

10° 

15° 

1.292 

I.290 

1.285 

1.276 

1.308 

1.308 

1.312 

1.314 

1.328 

1.329 

1.340 

1.349 

1.354 

1.357 

1.375 

• 1.391 

1.387 

1.391 

1.418 

1.443 

1.424 

1.433 

1.470 

I.505 

1.469 

1.481 

1.532 

1.582 

1.520 

1-537 

1.604 

1.672 

1.576 

1.598 

1.687 

1.780 

1.635 ' 

1.665 

1.782 

1.907 

•1.695 

1.736 

1.891 

2.058 

1.756 

1.810 

2.015 

2.241 

1.812 

1.886 

2.158 

2.469 

1.858 

1.957 ■ 

2.327 

2.766 

1.878 

2.015 

2.532 

3.182 

1.827 

2.018 

2.769 

3.802 

1.594 

1.836 

2.911 

4.692 

i.i4o 

1.306 

2.365 

4.891 

1.022 

1.000 

1.283 

2.514 

1.219 

1.133 

1.000 

1.189 

1.381 

. 1.294 

1.083 

1.000 

1.470 

1.396 

1.195 

1.055 

1.511 

1.449 

1.278 

1.137 

1.521 

1.471 

1.329 

1.204 

-1.513 

1.472 

1.358 

1.251 

1.493 

i.46o 

1.369 

1.281 

1.466 

i.44o 

1.368 

1.296 

1.436 

1.416 

1.359 . 

1.302 

l.4o6 

1.390 

1.346 

1.301 

1.377 

1.365 

1.332 

1.297 

1.351 

1.342 

1.318 

1.291 

1.330 . 

1.322 

1.306 

1.286 

1.313 

1.308 

1.297 

1.284 

1.301 

1.298 

1.294 

1.286 

1.296 

1.296 

1.296 

1.294 
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TABLE 6.- THE JOUKOWSEI AND ACTUAL EROF3LES 1 


[M M = 0.685; a ± = 0°] 



Actual: Z A = X A + iy A 
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Figure 1.- -Velocity distribution. M„ = 0.5; e = 0.15; R = 1.15. 
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Figure 2.- Velocity distribution. M«,= 0.685; e = 0.15; R = 1.15. 




